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ABSTRACT 

Using estimates of dark halo masses from satellite kinematics, weak gravitational lens- 
ing, and halo abundance matching, combined with the TuUy-Fisher and Faber- Jackson 
relations, we derive the mean relation between the optical, Vopt, and virial, V200J circu- 
lar velocities of early- and late-type galaxies at redshift z ~ 0. For late-type galaxies 
V^pt ~ V200 over the velocity range Vopt = 90 — 260 km s^^, and is consistent with 
Kjpt = VJnax.h (the maximum circular velocity of NFW dark matter haloes in the con- 
cordance ACDM cosmology). However, for early-type galaxies Vopt 7^ V200, with the 
exception of early-type galaxies with V^opt — 350 km s^^. This is inconsistent with early- 
type galaxies being, in general, globally isothermal. For low mass (Vopt ^ 250 km s""'^) 
early- types Vo-pt > Knax.h, indicating that baryons have modified the potential well, 
while high mass (T^pt ^ 400kms~^) early-types have Vopt < Knax,h- Folding in mea- 
surements of the black hole mass - velocity dispersion relation, our results imply that 
the supermassive black hole - halo mass relation has a logarithmic slope which varies 
from ~ 1.4 at halo masses of ~ lO^^/i~^M0 to ~ 0.65 at halo masses of lO^^-^ft.~^M0. 
The values of Vopt/V2m we infer for the Milky Way and M31 are lower than the 
values currently favored by direct observations and dynamical models. This offset is 
due to the fact that the Milky Way and M31 have higher Vopt and lower V200 com- 
pared to typical late-type galaxies of the same stellar masses. We show that current 
high resolution cosmological hydrodynamical simulations are unable to form galax- 
ies which simultaneously reproduce both the K)pt/V2oo ratio and the V^opt — -^'^star 
(TuUy-Fisher/Faber- Jackson) relation. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters 
- galaxies: haloes - galaxies: kinematics and dynamics ~ galaxies: spirals 



1 INTRODUCTION 

It is theoretically expected and observationally established 
that galaxies are surrounded by extended haloes of dark 
matter (White & Rees 1978; Blumenthal et al. 1984; van 
Albada & Sancisi 1986; Zaritsky & White 1994; Brainerd 
et al. 1996; Prada et al. 2003). In recent years much progress 
has been made in understanding the relation between the 
masses of dark matter haloes and the properties of the galax- 
ies that reside in them. The majority of work in the liter- 
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ature has focused on the relation between halo mass and 
galaxy luminosity or stellar mass (e.g. Yang et al. 2003; 
Kravtsov et al. 2004; van den Bosch et al. 2004; Hoekstra 
et al. 2005; Mandelbaum et al. 2006; Conroy et al. 2007; 
More et al. 2009). An alternative approach is to link galaxies 
to haloes via kinematics, for example, by measuring the rela- 
tion between the circular velocity (V'circ('') ~ \J GAI{< r)/r, 
for a spherical system) within the optical part (e.g. the half 
light radius or 2.2 disk scale lengths) of galaxies, U>pt, to the 
circular velocity at the virial radius of the dark matter halo, 

V200- 

This approach has the advantage that it is free from un- 
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certainties in luminosities or stellar masses due to uncertain- 
ties in extinction from dust, stellar populations or the stellar 
initial mass function (IMF). In addition, because it is a dy- 
namical link, it gives a direct measurement of the slope of the 
global mass density profile within the virial radius. For ex- 
ample, for a singular isothermal density profile, p(r) cx , 
and thus K = const, which implies Vopt/V2oo = 1- By mea- 
suring a departure from this prediction, one can rule out 
galax;ies being globally isothermal. This is a relevant issue 
because it is common practice to infer halo masses from 
observed Vopt by assuming Vopt = V200, even though this as- 
sumption lacks both theoretical and observational support. 
For example, in ACDM cosmologies dark matter haloes are 
not isothermal, so a detection of non-isothermality would be 
useful from a ACDM haloes perspective. 

The K)pt/V2oo ratio also contains information on the 
relative importance of baryons vs dark matter in the opti- 
cal regions of galaxies. For galaxy mass ACDM haloes, the 
ratio between the maximum circular velocity of the dark 
matter halo, Vmax.h, and the virial velocity of the dark mat- 
ter halo Vmax,h/V200 1-1 - 1-2 (Bullock et al. 2001), and 
thus if Vopt/V2oo is observed to be significantly larger than 
this, it indicates that baryons have modified the potential 
well. Using weak lensing measurements of halo masses com- 
bined with the TuUy-Fisher (1977, hereafter TF) relation, 
Seljak (2002) inferred Vopt/V2oo = 1.8 (with a 95% C.I. lower 
limit of 1.4), for late- type L* galaxies. Such high values of 
Vopt/V2oo are naturally explained by the combined effects of 
baryons adding to the optical circular velocity directly, and 
indirectly by inducing contraction of the dark matter halo 
(Blumenthal et al. 1986). 

However, a high value of V^pt/V2oo introduces problems 
in reconciling the TF relation with the halo mass function 
and galaxy luminosity function, LF. Semi-analytic galaxy 
formation models that are able to simultaneously reproduce 
the zero point of the TF relation and the galaxy LF as- 
sume that Vopt = V200 (e.g. Somerville & Primack 1999), 
or Vopt = Vinax,h (e.g. Croton et al. 2006). However, models 
that take into account the contribution of the baryons to the 
rotation curve and the effect of halo contraction have been 
unable to match both the LF and TF relation (e.g. Benson 
et al. 2003; Benson & Bower 2010). Dutton et al. (2007) 
showed that the Vopt/V2oo ratio can place strong constraints 
on dark halo structure. In particular they argued that mod- 
els with adiabatic contraction, standard ACDM halo con- 
centrations, and standard stellar IMFs are unable to si- 
multaneously reproduce the zero points of the TF and size- 
luminosity relations, and the requirement that Vo-pt — Vmax,h 
(in order to reproduce the LF). However, if halo contrac- 
tion was somehow avoided, or even if haloes could expand 
in response to galaxy formation, then models could be con- 
structed that reproduced all of the observational constraints. 

In this paper we derive the relation between Vo-pt and 
V200 for both early- (red, bulge dominated) and late-type 
(blue disk dominated) galaxies at redshifts 2; ~ 0. For 
late-types, we use Vopt = V2.2, where V2.2 is the rotation 
velocity at 2.2 disk scale lengths. For early-types we use 
Vopt = 1.65cr(7?5o), where a{R^o) is the velocity dispersion 
within the projected half light radius of the galaxy, and 
the factor of 1.65 is from Padmanabhan et al. (2004, see 
also ^2.4|l . We determine the mean relation between virial 
velocity and stellar mass using published measurements of 



halo masses from satellite kinematics, weak lensing, and halo 
abundance matching studies. By comparing these relations 
with the observed TF and Faber Jackson (1976, hereafter 
FJ) relations we derive the mean relation between Vopt and 
V200- In a future paper (Dutton et al. in prep) we will com- 
bine these results with other scaling relations of early- and 
late-type galaxies to place constraints on the structure of 
their dark matter haloes. 

This paper is organized as follows. In §2 we describe the 
observations of halo masses and galaxy scaling relations. In 
§3 we derive and discuss the relation between Vopt and V2oo- 
In §4 we discuss the implications of our results for the slope 
of the relation between super-massive black hole mass and 
halo mass. In §5 we compare our results with predictions 
for ACDM dark matter haloes. In §6 we compare our re- 
sults with halo masses derived for the Milky Way and M31 
galaxies, as well as galaxies formed in cosmological hydro- 
dynamical simulations. In §7 we give a summary. 

Unless otherwise specified, throughout this paper we 
adopt a Hubble parameter, Hq = 100 /ikms~^ Mpc~^, i.e. 
stellar masses are expressed in /i~'^Mq units, halo masses 
are expressed in /i~^M0 units, and halo sizes are expressed 
in /i~^kpc units. 



2 OBSERVATIONS 

This section gives an overview of the observational data 
we use to derive the relation between optical and virial 
velocities of galaxies. This paper primarily uses observa- 
tional data from studies based on the Sloan Digital Sky Sur- 
vey (SDSS), Data Release 4. The SDSS (York et al. 2000; 
Stoughton et al. 2002; Abazajian et al. 2004; Adelman- 
McCarthy et al. 2006) is an extensive photometric and spec- 
troscopic survey of the local universe. 

2.1 Virial Masses 

There are a number of techniques that are used to deter- 
mine virial massed of dark matter haloes. Direct observa- 
tional measurements of virial masses can be obtained with 
weak galaxy-galaxy lensing (e.g., Brainerd et al. 1996; Hud- 
son et al. 1998; Wilson et al. 2001; Guzik & Seljak 2002; 
Hoekstra, Yee & Gladders 2004; Mandelbaum et al. 2006; 
Cacciato et al. 2009) and satellite kinematics (e.g., Zaritsky 
& White 1994; Prada et al. 2003; van den Bosch et al. 2004; 
Conroy et al. 2005). A limitation of these techniques is that 
the signal for individual galaxies is usually too weak to give 
statistically significant measurements. This means that halo 
masses must be obtained by stacking many galaxies together 
(typically within some luminosity or stellar mass bin). An- 
other technique, known as the (sub) halo abundance match- 
ing method (Kravtsov et al. 2004; Vale & Ostriker 2004; 
Conroy et al. 2006; Conroy & Wechsler 2009), determines 
the relation between stellar mass and halo mass by assum- 
ing that there is a one-to-one mapping (sometimes scatter 
in this mapping is included) between the number density of 

^ Note that when we refer to the virial mass of a dark matter 
halo we implicitly include the baryonic mass as well as the dark 
matter mass 
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Table 1. Cosmological and stellar mass parameters adopted by studies of the halo masses - stellar mass relation. Notes — Col. (2) Method 
used to derive the halo mass - stellar mass relation: GC— Group Catalog; AM-Abundance Matching; WL-Weak Gravitational Lensing; 
SK-Satellite Kinematics. Col. (3) Mean density of dark matter haloes in units of the critical density of the universe (Avir = p/Pcrit)- 
Col. (4) Redshift zero cosmological matter density. Col. (5) Redshift zero cosmological dark energy density. Col. (6) Redshift zero Hubble 
parameter in units of kms~^ Mpc~^. Col. (7) Reference for stellar masses and stellar initial mass function: B03-Bell et al. (2003); 
BR07-Blanton & Roweis (2007); G06-Gallazzi et al. (2006); K03-Kauffmann et al. (2003); P07-Panter et al. (2007); diet-Sal.- Salpeter 
(1955) — 0.15dcx; Kroupa-Kroupa (2001); Chab.-Chabrier (2003); Col (8) Offset we apply to stellar mass, in dex. 
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dark matter haloes (including sub-haloes) and the number 
density of observed galaxies as a function of stellar mass. A 
related method, which we term the group catalog method, 
is to run a galaxy group finder and to assign each group 
a halo mass using abundance matching (this time not in- 
cluding sub-haloes) under the assumption that there is a 
one-to-one mapping between halo mass and the total stellar 
mass of all group members. 

In this paper we consider measurements based on all 
of these techniques: Weak lensing (WL) from Mandelbaum 
et al. (2006; 2008) and Schulz et al. (2010); Satellite kine- 
matics (SK) from Conroy et al. (2007); More et al. (2010); 
and Klypin, Prada & Montero-Dorta (2010, in prep); halo 
abundance matching (AM) from Moster et al. (2010), Guo 
et al. (2010) and Behroozi et al. (2010); Group catalogs 
from Yang et al. (2007). The SK measurements from Conroy 
et al. (2007) are somewhat different from what is reported 
in that paper. Here, we use stellar masses based on Bell 
et al. (2003), and we use halo masses derived from samples 
of SDSS galaxies that have not been randomly diluted in 
number density (as was done in order to mimic the DEEP2 
selection function). Otherwise the sample and methodology 
is identical to Conroy et al. (2007). 

2.1.1 Homogenization 

All of these methods are cosmology dependent to some ex- 
tent. Deriving halo masses from weak galaxy-galaxy lensing 
or satellite kinematics requires some knowledge of the struc- 
ture of the dark matter halo. The halo abundance matching 
method requires a halo mass function (which depends di- 
rectly on a number of cosmological parameters). 

The virial masses were calculated under different as- 
sumptions and cosmologies. See Table [T] for the values of 
virial radius overdensity (Avir), redshift zero matter density 
(f2m), redshift zero dark energy density (Oa), redshift zero 
Hubble parameter {Ho). Note that for Yang et al. (2007) 
the halo masses are based on the halo mass function of 
Warren et al. (2006). The conversion between this mass 
function and spherical overdensity masses is non-trivial, as 



the Avir is dependent on resolution. Here we assume that 
the Warren et al. (2006) halo mass function corresponds to 
Avir = SOOpmcan- While it is not possible to bring the mea- 
surements from each author onto exactly the same cosmol- 
ogy, we can at least adopt the same definition of halo mass. 
For the results in this paper, the choice of halo definition is 
arbitrary, so we adopt the most convenient definition. 

Halo masses are commonly defined via Avir, the over- 
density of the halo with respect to the critical density of the 
Universe: 

{p) = Mvir/^7ri?^ir = Avir Pcrit ■ (1) 

Here Mvir is the virial mass, i?vir is the virial radius, and 
Pcrit ~ 3Hq/8tvG is the critical density of the Universe. 

In this paper we adopt Avir ~ 200, so that halo masses 
are given by M200, halo sizes by 7i2oo, and halo virial veloc- 
ities by V200. This definition results in a simple conversion 
between halo virial size, virial velocity, and virial mass at 
redshift 2 = 0: 

, V200 , -R200 1 , f ^ M200 \ ,„\ 
= [ft^ = 3 (^^ [fe^TM^ ) ' (2) 

where G ~ 4.301 x 10"*^ ( kms-^)^ kpc Mq-^ This definition 
also has the advantage that halo masses are independent of 
the matter density of the universe (which is still subject to 
a significant uncertainty). 

The conversion between the halo masses defined with 
different values of Avir depends on the concentration of the 
halo. For the halo masses from Mandelbaum et al. (2008) 
and Schulz et al. (2010) we use the halo concentrations that 
these authors derive from the weak lensing fits. For the other 
authors we adopt the halo concentration - mass relation for 
relaxed haloes in a WMAP 5th year cosmology (WMAP5, 
Dunkley et al. 2009), from Maccio et al. (2008): logio c = 
0.830 - 0.098 logi(j(M2oo/lO^^/i-^M0). 

Stellar masses were calculated using different IMFs and 
assumptions about the star formation histories and metal- 
icities. While we cannot account for differences in modeling 
methods, we can attempt to convert stellar masses to a uni- 
form IMF. As our fiducial stellar masses we adopt those from 
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Figure 1. Relation between average virial mass at fixed stellar mass, (M20o)i a-iid stellar mass, Mstar, expressed in terms of (Af2oo)/Afatar 
vs Mstar, using data from the literature as indicated. The various methods (satellite kinematics, SK, weak gravitational lensing, WL, halo 
abundance matching, AM, and group catalogs, GC) all yield consistent results within the quoted uncertainties. All error bars are 95% or 
2a. The black lines and points show the relations derived from taking all types of galaxies and are shown in both panels. The red points 
show the relation derived for early- type galaxies (left panel), and the blue points show the relation for late- type galaxies (right panel). 
The solid lines show our fitting formula (determined by eye) for the mean and 2a uncertainty in the relation between (M200} (Mstar) and 
Mstar as given by Eq. [S] The numbers on the right hand y-axis of each panel indicate the percentage of the cosmological baryons that 
reside in stellar mass in central galaxies (assuming h = 0.7 and /bar = f^b/f^m = 0.165). 



Bell et al. (2003) minus 0.1 dex. Unless otherwise stated, we 
assume that these masses correspond to a Chabrier (2003) 
IMF. Our adopted conversions between the stellar masses 
calculated by various authors to this definition are given in 
Table [T] We note that the dilTerences between various stel- 
lar mass estimators vary with stellar mass, even when the 
same IMF is adopted. For simplicity we have chosen con- 
versions most appropriate for stellar masses in the range 
Mstar =^ 10^° - 10" /i"^M0. The differences between masses 
from Bell et al. (2003), Kauffmann et al. (2003), and Blan- 
ton & Rowels (2007) are based on the comparison made in 
Appendix A in Li & White (2009). 

For the Mandelbaum et al. (2008) results the halo 
masses are given for r-band luminosity bins. We convert 
these luminosities into stellar masses by using the Faber- 
Jackson relations from Gallazzi et al. (2006), which yield 
the following mean relation between stellar mass and r-band 
luminosity: logjg Mstar = 1.093 log j^q I/r — 0.573. 



2.2 Relation Between Halo Mass and Stellar Mass 

Fig. [T] shows the relation between halo virial mass, M200, 
and galaxy stellar mass, Mstar, as expressed in terms of 
(M2oo)(Mstar)/Mstar VS Afstar, for the different methods de- 
scribed in i ]2.1l There are important differences between 
the way the relation between halo mass and stellar mass 
is usually expressed between the various methods. The AM 
method usually gives the mean of the log of stellar mass 
as a function of halo mass: (logjQ A/star) (Af20o)- While the 



WL and SK masses give the log of the mean halo mass as 
a function of stellar mass: logj^Q{A/2oo)(Mstar)- For low stel- 
lar masses these measurements give equivalent relations be- 
tween halo mass and stellar mass. However, at high masses, 
these two definitions give diverging relations between halo 
mass and stellar mass, with log]^Q(M2oo)(Mstar) giving sys- 
tematically lower halo masses than (logjQ Mstar) (M200) (e.g. 
see Fig. 10, Behroozi et al. 2010). This effect is caused by a 
combination of the following: scatter in stellar mass at fixed 
halo mass; the number density of haloes strongly decreases 
with increasing halo mass; and the relation between stellar 
mass and halo mass is shallow at high masses (e.g. More, 
van den Bosch & Cacciato 2009). Since the primary goal 
of this paper is to determine the 14pt/V2oo ratio as a func- 
tion of observables such as Vopt and Mstar, we thus focus on 
logio(M2oo) (Mstar). For the AM masses (Moster et al. 2010; 
Guo et al. 2010; Behroozi et al. 2010) we use (M200) (Mstar) 
provided to us by the authors. These have been derived as- 
suming a scatter of 0.15, 0.20, and 0.15 dex in stellar mass 
at fixed halo mass. For the GC masses (Yang et al. 2007) 
we use the (M200) — Mstar relation as calculated by More 
et al. (2010). 

There are different ways of splitting galaxies into the 
broad categories of early- and late-type. Some studies use 
the bimodality seen in the color-magnitude diagram (with 
early-types being red, late-types being blue), while others 
divide based on the concentration of the light profile (with 
early-types being high concentration, late-types being low 
concentration). In this work we assume the two methods 



© 2010 RAS, MNRAS OOO.ITlfTTl 



Connecting Galaxies with Dark Matter Haloes 5 



Table 2. Parameters of double power-law fitting formula (Eq. [Sjl 
to the y = (M2oo>/Afs tar vs X — A/star relations in Fig.^ Masses 
are in units of h = 1. 
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are equivalent, although in detail they are unlikely to be so. 
For example, red galaxies can be contaminated with dusty 
spirals (which are intrinsically blue) , and high concentration 
galaxies can be contaminated with blue spiral galaxies with 
significant bulges. The black lines and points in Fig. [T] show 
the measurements derived from taking all types of galax- 
ies (both panels). The red points show measurements for 
early-type galaxies (left panel), while the blue points show 
measurements for late- type galaxies (right panel). 

Since late-type galaxies dominate at low masses, and 
early-type galaxies dominate at high masses, the relation 
between A/2ao and A/star for early-type (late-type) galax- 
ies should converge to those using all galaxy types at high 
(low) stellar masses. While this is the case at high masses, 
at low masses the AM method gives halo masses that are 
systematically lower, by a factor of ~ 2, compared to the 
SK method. We attribute this discrepancy to an unknown 
systematic uncertainty in either one of or both the AM or 
SK methods. 

We use the following function for y{x) to provide a fit- 
ting formula for the y = (M200) /Avatar vs. x — Afatar relation 
for early- and late-type galaxies: 
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Here a is the logarithmic slope at a; <C so, /? is the logarith- 
mic slope at X ^ Xo, Xo is the transition scale, yo — y{xo) 
is the y value at the transition scale, and 7 controls the 
sharpness of the transition. 

Given the diverse range of measurement techniques and 
error estimation, we fit the relation between {Ai^20o)/Afatar 
and Afstar by eye. The parameters of our fits are given in 
Table (2] The upper and lower limits were chosen to bracket 
the error bars of the measurements (which are 2a or 95% 
C.I.), as well as taking into account the halo abundance 
matching results for low mass late-types. Thus these limits 
can be thought of as roughly 2g systematic errors. Note that 
we do not include the WL results for low mass early types 
from Mandelbaum et al. (2006) in our fits, as these have 
been superseded with more recent measurements (Mandel- 
baum et al. 2008; Schulz et al. 2010). These more recent 
papers use isolation criteria to select central galaxies, rather 
than selecting all galaxies (of a given type) and then trying 
to model the weak lensing signal in terms of centrals and 
satellites, as was done in Mandelbaum et al. (2006). 



2.2.1 Galaxy Formation Efficiency 

The numbers on the right vertical axes in Fig. [1] show the 
percentage of the cosmologically available baryons that end 
up as stars in a given galaxy, where we adopt a cosmological 
baryon fraction of /bar = flb/f^m = 0.165, and a Hubble 
parameter oi h = 0.7. We refer to this as the integrated 
star formation efficiencjO, esF ~ Afstar/(/barA^20o). When 
including the cold gas, this parameter becomes the galaxy 
formation efficiency ecF = (Afstar + Afgas)/(/barAf20o). For 
early-type galaxies we assume that the cold gas fractions are 
small enough to be ignored, and thus ecp = esF- 

For late-type galaxies the integrated star formation 
efficiency increases with increasing stellar mass. Ranging 
from esF = 8.9^5^1"% (2cr) at a stellar mass of A/star = 
lO^-^/i'^M©, to esF = 2Q.2,t\ol% (2o") ^t a stellar mass 
of A/star ~ lO^^/i~^M0. Using the relation between cold 
gas fraction (both atomic and molecular) and stellar mass 
from Button & van den Bosch (2009): /gaa = 0.374 - 
0.162(log]^Q Afstar — 10), the corresponding galaxy formation 
efficiencies are ecF = 16.7tft^% (2cr)> ^gf = 2,?,.?.trA% 
(2a). 

For early-type galaxies the galaxy formation efficiency 
peaks at ecF = 12.4l5°g^% (2cr), at a stellar mass of 
Afstar = lO^° ''/i"^M0. Below this mass the galaxy forma- 
tion efficiency is consistent with being constant. Above this 
mass the galaxy formation efficiency decreases, such that at 
the highest stellar masses probed, Afstar = 1O""/i-2M0, 
esF — 2.8l'^ g% {2a). Note that this galaxy formation ef- 
ficiency is only counting the central galaxy. These massive 
early-types are centrals in clusters, and so there is a sub- 
stantial (possibly dominant) contribution to the total stellar 
mass budget in those haloes from satellite galaxies (e.g. Lin 
& Mohr 2004). 

Using weak galaxy-galaxy lensing derived halo masses 
Hoekstra et al. (2005) found an average galaxy formation 
efficiency of ~ 33% for blue galaxies (taking into account 
the cold gas fraction) and 14% for red galaxies. These results 
are consistent with those presented here. However, we note 
that since the galaxy formation efficiencies vary with stellar 
mass, the average value will depend on the sample selection, 
and it is therefore not very meaningful to quote a single 
value. 

2.3 The Tully-Fisher Relation 

The original Tully-Fisher relation was between B-band lumi- 
nosity and HI linewidth. The term TF relation is currently 
used to describe the broader class of relations between lumi- 
nosity (at optical to near-IR wavelengths) or mass (stellar 
or baryonic) and rotation velocity (measured in a number 
of different ways) for late-type galaxies. As such there is no 
universal, or correct, definition of the TF relation. Often 
different versions are better suited to different applications. 

Since the measurements of halo mass that we use in 
this paper are as a function of stellar mass, here we focus 
on the stellar mass TF relation. For rotation velocities our 

^ Note that what we term the integrated star formation efficiency 
should not be confused with the integral of the star formation 
history of a galaxy, which will be larger than esF due to mass lost 
through stellar winds and supernovae. 
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Figure 2. Stellar mass - TuUy-Fisher relations derived from the 
Pizagno et al. (2007) and Verheijen (2001) data. The thick lines 
show fits of the form y = a + b{x — 10), with best fit a and b as 
indicated. The TF relations from Bell & de Jong (2001), Pizagno 
et al. (2005) and Button et al. (2007) are given by the dotted- 
short-dashed, long dashed, and short dashed lines, respectively. 
For comparison purposes, the dotted line shows the Faber- Jackson 
relation from Gallazzi et al. (2006), after correcting velocity dis- 
persions to circular velocities assuming VcircC^so) = 1.65(T(iJ5o). 

favoured definition is V2.2, the rotation velocity measured at 
2.2 I-band exponential disk scale lengths. We adopt this for 
several reasons: The main reason is that the velocity is mea- 
sured at a well defined radius, which enables mass models 
to be fitted to the TF relation; Secondly, this radius is small 
enough that baryons and dark matter contribute roughly 
equally to the circular velocity, and thus V2.2 provides a 
velocity that is sensitive to the distribution of baryons and 
dark matter in galaxies. Lastly, V2.2 is relatively easy to mea- 
sure observationally from low to high redshifts, and thus the 
V2.2 — Afatar relation can be measured for large statistically 
complete samples of galaxies. 

Other velocity measures, such as Vmax, the maximum 
observed rotation velocity, or Vflat, the rotation velocity in 
the flat part of the rotation curve, may appear simpler or 
more fundamental, but since there is no radius associated 
with the rotation velocity, these definitions are not as use- 
ful for mass models. Furthermore, since Vnax is measured at 
large radius, it is insensitive to the distribution of baryons 
and dark matter on small scales, and is observationally chal- 
lenging to measure. 

In this paper we base our stellar mass TF relation on the 
sample of ~ 160 SDSS galaxies from Pizagno et al. (2007; 
hereafter P07). These galaxies were not selected to be late- 
type, although the requirement that they have sufficiently 
strong and extended Ha emission to yield a useful rotation 
curve is equivalent to a late- type selection based on color. 
We compute stellar masses for the P07 galaxies using the 
r-band mass-to-light ratio vs g — r color relation from Bell 
et al. (2003), with a -0.10 dex offset (corresponding to a 



Chabrier IMF): log^o Afstar/ir = -0.506 + 1.097(ff - r). In 
order to minimize inclination uncertainties and extinction 
corrections we restrict our sample to galaxies with axis ratios 
0.3 < b/a < 0.5. These data are shown as black filled circles 
in Fig. [2I The stellar mass TF relation that we derive here 
from the P07 data is shown as the black line in Fig. [5] and 
is given by: 

1°^-' Ik^ = 2-143 + 0.281 (log,, p,^^) . (4) 

The uncertainty on the zero point is 0.027, the uncertainty 
on the slope is 0.003, and the intrinsic scatter is 0.05 dex in 
V2.2- For comparison, the V2.2 — Mstm TF relation derived 
in Button et al. (2007), using the much larger but more 
heterogeneous sample from Courteau et al. (2007), is shown 
as the black short dashed line, and given by: 

Ik^ = + [10^^-^^Mel ) • 

The uncertainty on the zero point is 0.039, the uncertainty 
on the slope is 0.004. This relation has a shallower slope, 
but Fig. [5] shows that the differences between these two TF 
relations are very small. Fig. [2] also shows the V2.2 — Afstar 
TF relation from Pizagno et al. (2005) which was based on 
a subset of 81 disk dominated galaxies. This relation has a 
similar slope and zero-point to the one derived here. 

The red open circles show the Vmax — Mstar relation 
derived from data in Verheijen (2001). Stellar masses are 
calculated using the relation between 7-band stellar mass- 
to-light ratio and {B — R) color from Bell et al. (2003) , with 
a 0.1 dex offset (corresponding to a Chabrier (2003) IMF): 
logio Matar/L/ = -0.505 0.518(5-7?). We adopt the HST 
Key Project distance to the Ursa Major Cluster of 20.7 Mpc 
(Sakai et al. 2000), which was also used by Bell & de Jong 
(2001), but larger than the distance of 18.6 Mpc used by 
Verheijen (2001). To convert masses to units we adopt 
h = 0.7. A fit to the data is shown by the red dot-long- 
dashed line and is given by: 

11^ = 2-185 + 0.235 (log p,^^). (6) 

The uncertainty on the zero point is 0.057, the uncertainty 
on the slope is 0.006. This relation has a shallower slope 
than the V2.2 — Mstar relation. At high stellar masses Knax — 
V2.2, and thus the difference is driven by Knax being higher 
than V2.2 for low mass galaxies, which indicates that rotation 
curves are still rising at 2.2 disk scale lengths in low mass 
(Vmax — lOOkms"^) galaxies. 

For completeness we also show (dot-short-dashed line 
in Fig. 01 the Vfo^t - Mstar TF relation from Bell & de Jong 
(2001), which also used data from Verheijen (2001). We show 
the relation using stellar masses derived from /-band lumi- 
nosities, {B — R) colors, mass dependent extinction correc- 
tions, and h = 0.7. This relation has a similar slope and zero 
point to the Vmax — Mstar rclatiou that we derive from the 
same data set. 

2.4 The Faber-Jackson Relation 

Gallazzi et al. (2006) measured the FJ relation for early- 
type galaxies using SDSS data. Velocity dispersions were 
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corrected from the SDSS 3.0 arcsec diameter aperture to 
the projected r-band half hght radius, R50, by the standard 
assumption that the radial profile of the velocity dispersion 
has a log slope of —0.04 (j0rgensen 1999). Note that these 
corrections are typically of order a percent. Stellar masses 
(from Gallazzi et al. 2005) were computed using a Chabrier 
IMF, and assuming h = 0.7. The stellar mass FJ relation is 
given by: 

The relation is valid for stellar masses in the range 10^ < 
Mstar < 10" ''/i"^Mq . The uncertainty on the slope is 0.020, 
and the scatter is 0.071 dex in cr(i?5o). 

We adopt the conversion factor between velocity disper- 
sion within the projected half light radius, cr(7?5o), and circu- 
lar velocity (at the projected half light radius) as derived in 
Padmanabhan et al. (2004): V50 = Khc(-R5o) = 1.65 cr(i?5o). 
This has an uncertainty of about 10% depending on the 
anisotropy profile. For comparison, Seljak (2002) assumed 
Kirc(i?5o) = l.5a{R5o), whereas Wolf et al. (2009), argue 
that = •^'^iob'''i/2/G, independent of the anisotropy. 

Here Mi/ 2 is the mass enclosed within a sphere of radius 
ri/2, the deprojected half light radius, and aios is the line 
of sight velocity dispersion of the system. Recasting this in 
terms of circular velocity: V1/2 = \/3o"ios — 1.73(Jios, which 
is within 5% of our adopted conversion. 

The principle motivation for this conversion is to put ro- 
tation velocities and velocity dispersions on the same scale, 
so that the ratio between the optical and virial velocities 
of early- and late-type galaxies can be more directly com- 
pared. The relation between V50 and stellar mass for early 
types is shown as a dotted line in Fig. [J] This relation has 
a very similar slope to that of the V2.2 — Mstar TF relation 
for late-types, but has a higher normalization (in velocity) 
by ~ 0.13 dex. 

Since we are using different definitions for the optical 
circular velocity for early and late-types it is worth dis- 
cussing how fair it is to compare them directly. Recall that 
for late-type galaxies we use circular velocity measured at 
2.2 I-band disk scale lengths, whereas for early-type galaxies 
we use circular velocity measured at the projected half light 
radius. For a bulge-less exponential disk, 65% of the light is 
enclosed within 2.2 disk scale lengths. For a bulge plus ex- 
ponential disk, the fraction of enclosed light will typically be 
higher. For an exponential disk the maximum circular veloc- 
ity occurs at 2.16 disk scale lengths, or 1.29 half mass radii. 
For a Hernquist sphere (which is close to a deprojected de- 
Vaucouluers sphere) the maximum circular velocity occurs 
at 0.551 projected half mass radii, or 0.414 3D half mass 
radii. 

Thus relative to the half light radius we are mea- 
suring Vopt at a slightly larger radius for late-types than 
early-types, but relative to the peak circular velocity of the 
baryons, we are measuring T4pt at a smaller radius in late- 
types than early-types. However, since both early and late 
type galaxies have roughly constant circular velocity pro- 
files in the optical parts of galaxies (e.g. Rubin et al. 1985; 
Koopmans et al. 2006), the small differences between our ve- 
locity definitions for early and late-types does not introduce 
a significant bias in the Vopt — A/star relations. 



The different normalizations of the TF and FJ rela- 
tions has important consequences for the universality of the 
M200 — Mstar and V200 — Vopt relations. It implies that at 
most one of these two relations can be universal (i.e. true 
for all types of galaxies). Current observations suggest that 
both the M200 - A'/star (Fig. [T]) and V200 - Kpt (see ^Hl be- 
low) relations are different for early- and late-types. How- 
ever, given the systematic uncertainties on halo masses, we 
cannot rule out the possibility that either of these relations 
is universal. A more definitive answer will require more ac- 
curate measurements of halo masses as a function of galaxy 
type. 



3 RELATION BETWEEN OPTICAL AND 
VIRIAL CIRCULAR VELOCITIES 

Given a relation between virial mass and stellar mass, we 
can trivially compute a virial velocity - stellar mass relation 
using Eq. (2] Using the fitting formula to the data in Fig. [T] 
these relations for early- and late-type galaxies are shown 
in the upper panels of Fig. [31 By comparing these relations, 
with the FJ (Eq. (T]) and TF (Eq. ^ relations (dashed lines 
in Fig. [3]), we can derive the average ratio between optical 
circular velocity (V50 for early-types, and V2.2 for late-types) 
and virial circular velocity (lower panels in Fig. [3]). Note 
that the dominant uncertainty in the V^pt/Vjoo ratio is the 
uncertainty in V2oo. Uncertainty in the conversion from a to 
Vopt, which is just ~ 0.04dex (black shaded region in upper 
left panel of Fig. [3]), contributes negligibly, and is therefore 
ignored in our error analysis. 

For this procedure to be valid requires that the galaxies 
in the different data samples (e.g. TF vs M200 /Mstar) are 
representative of the same population, and that the stellar 
masses are comparable. We have verified that we get consis- 
tent results for K)pt/V2oo vs. Vopt when using the velocity - 
luminosity relations (i.e. V200 - Lr and Vopt - Lr) instead of 
velocity - stellar mass relations. Thus the main systematic 
uncertainty in our derivation of Vopt/V2oo is possible sample 
selection differences between TF and halo mass measure- 
ments. 

For late-type galaxies the virial circular velocity - stel- 
lar mass relation (V200 vs Afstar) is shown as blue shaded 
region in the top right panel of Fig. |3] There is some curva- 
ture to the V200 — Mstar relation, but within the uncertain- 
ties the data is consistent with a power-law over the stellar 
mass range Afstar = 10^" - lO" °/i"^M0. The slope and zero 
point of this relation closely resembles the optical velocity - 
stellar mass relation (dashed lines, Eq|4ll. This implies that 
Vopt/V2oo — 1, as is shown in the lower right panel. 

For early-type galaxies the V200 — Afstar relation is not 
well described by a single power-law (upper left panel in 
Fig. |3}. At low stellar masses the slope is similar to that 
of the FJ relation, but at high stellar masses the slope is 
much steeper than that of the FJ relation. This implies that 
V50/V200 is a strong function of stellar mass (lower left panel 
of Fig. E}. At a stellar mass of Mstar ~ 10"/i"^Mq, V50 ^ 
V200. For higher stellar masses, V50 < V200, while for lower 
stellar masses V50 > V200. 
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Figure 3. Circular velocity vs stellar mass relations (upper panels) for early- (left panels) and late-types (right panels). The shaded 
regions correspond to 2cr uncertainties. The stellar mass TuUy-Fisher (TF) and Faber-Jackson (FJ) relations are given by the dashed 
lines. For the FJ relation we have converted velocity dispersions into circular velocities assuming V50 = 1.65 (7(-R5o). The uncertainty in 
this conversion is given by the black shaded region, and is small compared to the uncertainty in halo virial velocities. The lower panels 
show the difference between the TF and FJ relations and the halo virial velocity relations. The points with error bars (2(t) show the 
values of V'opt/V200 derived for L* galaxies by Seljak (2002). 



3.1 Comparison with the Literature 

Seljak (2002) used the weak lensing virial mass-to-light ra- 
tios from Guzik & Seljak (2002) and the same method as 
used here to derive a l/opt/1^200 ratio for L* galaxies. For 
late-type L* galaxies (Kpt = 207kms"^) Seljak (2002) 
found Vopt/V2oo = 1-8 with a 2ct lower limit of 1.4, there 
was no upper limit due to the lensing masses being consis- 
tent with zero. For early-type L* galaxies [a — 177 km s^^) 
Seljak (2002) found V'opt/F200 = 1.68 ± 0.2 (2a uncer- 
tainty), assuming Vopt = 1.5(t. Using our adopted value of 
Vopt = 1.65(7 the result from Seljak (2002) corresponds to 
Kjpt/y200 = 1.85±0.22. These results are shown as filled cir- 
cles and error bars in Fig. [S] For both early- and late-types 
the results from Seljak (2002) are higher than our values, 
with the differences being larger than 2a. We trace this dis- 
crepancy to the relatively low values of M200/L used by 
Seljak (2002). These values are only marginally consistent 



with weak lensing results of Mandelbaum et al. (2006) which 
used a much larger sample than Guzik & Seljak (2002), as 
well as improvements in the modelling and calibration of the 
lensing signal. Furthermore, the results from Mandelbaum 
et al. (2006) have also been superseded with results from 
Mandelbaum et al. (2008) and Schulz et al. (2010) which 
find even higher halo masses. 

Seljak (2002) used the high value of Vopt/V2oo as evi- 
dence for adiabatic contraction of the dark matter halo. In 
light of the significantly lower value of Vo-pt /V200 that we find 
here, the conclusions of Seljak (2002) need to be revisited, 
which we do so in a future paper (Dutton et al. in prep). 

Eke et al. (2006) used the group abundance matching 
method with the 2dF Galaxy Redshift Survey Percolation- 
Inferred Galaxy Group (2PIGG) catalogue to derive a re- 
lation between halo mass (or equivalently halo circular ve- 
locity) and total &,7-band luminosity. They then convert to- 
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tal group luminosity into average central galaxy luminosity. 
By comparing their halo circular velocity - central galaxy 
luminosity relation with a B-band TF relation (from Bell 
& de Jong 2001) they derive a relation between Vflat/K-ir 
as a function of &j-band luminosity. Here Vkat is the ro- 
tation velocity in the flat part of the rotation curve, and 
K-ir is the circular velocity at the virial radius, defined with 
Avir — 100 (see Eq. [l}. For 6j-band luminosities in the 
range L^j = 10^'^ - lO^-^h'^LQ, Eke et al. (2006) found 
Vfiat/Vvir — 0.9 — 1.0. While for more luminous galaxies 
Vkat/K-ir declines with increasing luminosity. 

A direct comparison between our results and those 
of Eke et al. (2006) is complicated by the fact that Eke 
et al. (2006) derived Vkat/Vvir by comparing a Vvir — L 
relation for all types of galaxies with the TF relation for 
spiral (late-type) galaxies. There are also differences be- 
tween velocity definitions, but these are small (~ 5%) and 
mostly cancel out. With these caveats in mind, the values 
of Vflat/K-ir found by Eke et al. (2006) are consistent, al- 
beit on the low side, with our results for late-type galaxies. 
The lower values of Vflat/Vvir found by Eke et al. (2006) 
can be mostly explained by the fact that they compared 
an extinction un-corrected Vvir — Lb relation to an extinc- 
tion corrected Vaat — Lb relation. This inconsistency biases 
Vflat/Kir low by ~ 0.07 dex. 

Guo et al. (2010) use the relation between halo mass 
and stellar mass from abundance matching to construct a 
Knax.h— Mstar relation. They compare this to the Vflat— Mstar 
relation of spiral galaxies from Bell & de Jong (2001) find- 
ing that in general Vaat Knax.h- Thus their implied values 
of Vflat /V200 are slightly higher than we find here, but they 
are within the uncertainties on halo masses. Since we have 
also used a TF relation from Bell & de Jong (2001), the dif- 
ferences between our results and theirs must originate from 
differences in the halo mass - stellar mass relation. These dif- 
ferences can be seen in Fig.[T] Below a stellar mass of Mstar = 
lO^°/i"^M0 the Guo et al. (2010) relation is offset to lower 
halo masses, which explains the higher Vflat /V200 implied by 
their results. At a stellar mass of Mstar lQl^°-^h-'^MQ (cor- 
responding to M200 ^ lO^^ °/i"^M0 and V200 = 163 km s"^) 
the relations cross over. For galaxies with Vopt — 220kms~^ 
Guo et al. (2010) should have found V^pt < Vmax.h- Their 
result that V^pt ^ V^max.h was biased by the fact that the 
observed TF relation data that they adopted had only three 
galaxies with Vopt > 200kms~^, which happened to be bi- 
ased to high V^pt at fixed stellar mass (see Fig. [2]). 

Guo et al. (2010) claim that the reasonable agreement 
between the observed Vflat — Mstar TF relation and their 
derived V^max.h — Mstar TF relation implies that the ACDM 
cosmology does seem able to reproduce observed luminos- 
ity functions and TuUy-Fisher relations simultaneously. As 
shown below in !j5]we agree with the conclusion that V^pt — 
Vmax.h (for galaxics in haloes with V200 — 100 — 300 kms~^). 
However, the fact that one can simultaneously reproduce 
the LF and TF relations in ACDM if one sets Vopt = Vmax,h 
has been known from previous studies (e.g. Cole et al. 2000; 
Benson et al. 2003; Croton et al. 2006). The problem is that 
Vopt — Vmax.h is uot uaturally reproduced by cosmologi- 
cal simulations of disk galaxy formation (see ij6] below) or 
analytic models (Button et al. 2007). Rather than being a 
fundamental challenge for the ACDM paradigm, the most 
likely explanation for this difficulty is that we still lack an 



adequate understanding of galaxy formation, and inparticu- 
lar the response of dark matter haloes to galaxy formation. 



3.2 Are Early- Type Galaxy Mass Density Profiles 
Isothermal? 

It is well known that late-type galaxies often have roughly 
flat rotation curves over the radii probed by gas kinemat- 
ics. More recently Koopmans et al. (2006, 2009) used a joint 
lensing and dynamics analysis of a few dozen strong gravita- 
tional lenses from the SLAGS survey (Bolton et al. 2006) to 
measure the slope of the total (i.e. dark and baryonic) den- 
sity profile within the half light radius of early-type galaxies. 
These authors found that the total density slope was very 
close to isothermal. Furthermore, by combining the results 
from the SLAGS lenses with those from the LSD survey 
(Treu & Koopmans 2002), Koopmans et al. (2006) find that 
the total density slope is close to isothermal out to a redshift 

2 = 1. 

The range of radial scales probed by these studies is 
small: from 0.2 to 1.3 galaxy half light radii. It is thus possi- 
ble that the galaxies studied by Koopmans et al. were only 
locally isothermal over the small range of radii that they 
could probe. To obtain a measurement of the total mass den- 
sity slope over a larger range of radii, Gavazzi et al. (2007) 
measured the weak lensing signal from 18 strong lenses from 
SLAGS, and found that the total mass density profile is 
roughly isothermal out to large (few 100 kpc) radii. 

The lower left panel of [3] shows that while constant 
circular velocity profiles cannot be ruled out (at the 2a 
level) the majority of early-type galaxies favor V^pt 7^ V200- 
However, at a circular velocity of V^pt — 350kms~^ (i.e. 
a ~ 215kms"\ or Mstar ^ lO^/i'^M©) the Vopt/V2oo ratio 
is unity, which is a necessary, but not sufficient, condition 
for a galaxy mass distribution to be globally isothermal. The 
lenses used in Gavazzi et al. (2007) had velocity dispersions 
200 < cr < 330 km s"^ (i.e. log^o Vopt =i 2.6 ± 0.1). Thus our 
results are consistent with those of Gavazzi et al. (2007), but 
suggest that early-type galaxies can only be globally isother- 
mal (i.e. a constant circular velocity, V^irc, from the optical 
half light radius, -R50, to the virial radius of the host dark 
matter halo, R200), if at all, over a narrow range of velocity 
dispersions. Further support for the conclusion that early- 
type galaxies are not globally isothermal comes from Klypin 
& Prada (2009) who studied the kinematics of satellites, at 
distances of 50-500 kpc, around isolated red galaxies in the 
SDSS. 

We remind the reader that, while WL and SK analyses 
assume NEW profiles (which are not isothermal) to derive 
halo masses, the resulting halo mass is quite insensitive to 
the actual shape of the halo mass profile (e.g. van den Bosch 
et al. 2004; Gonroy et al. 2007). This is also supported by 
the fact that halo masses from WL and SK analyses are 
consistent with those from halo abundance matching, which 
makes no assumption about the structure of dark matter 
haloes. Furthermore, measurements of Vopt are independent 
from those of V200- This means, for example, that even if we 
assumed that haloes are isothermal in deriving halo masses 
from SK and WL, this would not guarantee that Vopt / V200 = 
1. Thus, there is therefore no circularity in our conclusion 
that the mass profiles of early-type galaxies, from the optical 
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half light radius to the virial radius of the halo, are in general 
non-isothermal. 



4 IMPLICATIONS FOR THE BLACK HOLE 
MASS - DARK HALO MASS RELATION 

There is a well studied correlation between the masses of 
central super-massive black holes (BH) and the bulge mass 
or velocity dispersion: the Mbh — o relation (Magorrian 
et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
Tremaine et al. 2002). There have been attempts to extend 
this relation to one between black hole mass and dark halo 
mass (e.g. Ferrarese 2002; Bandara, Crampton, & Simard 
2009), by or Vopt — V200- However, 

since this conversion is at best applicable as an average for 
samples of galaxies, it is currently not possible to determine 
the scatter in the black mass - halo mass relation. Here we 
focus on the implications of our results for the slope of the 
relation between black hole mass and halo mass. 

By using our relation between halo mass and stellar 
mass for early- type galaxies (Eq. |3] Table [2]), combined with 
the FJ relation (Eq. (Tjl we can derive a relation between 
halo mass and galaxy velocity dispersion. By combining this 
relation with the A/bh — o relation for elliptical galaxies 
from Giiltekin et al. (2009): logn, Mbh = (8.23 ± 0.08) + 
(3.96 ± 0.42)logio(cr/2QQkms"^), valid for 3 x IQ'^Mq < 
Mbh < 3 X IQ^Mq, we derive a relation between the mean 
halo mass as a function of black hole mass. This relation is 
shown in Fig. 31 where the shaded region corresponds to the 
2a uncertainty in the M200 — Mstar relation. The parameters 
of our fitting function are given in Tabled 

Observational and theoretical studies of the Mbh — M200 
relation often make the simplifying assumption that the ob- 
served circular velocity (in the optical part of a galaxy) is 
proportional to the circular velocity at the virial radius: 
V^jpt = 7V200, where a singular isothermal sphere corre- 
sponds to 7 = 1. For example Croton (2009) use this as- 
sumption, with 7 = 1, to derive a relation between black 
hole mass and halo mass Mbh oc Ml^f^^'^'^. Bandara 
et al. (2009) use strong lensing masses from the SLAGS 
survey (Bolton et al. 2006), and the implicit assumption of 
7 = 1, to derive a relation between black hole mass and halo 
mass, A^BH oc MIi^q^°-^^ . Adopting V200 = Kpt = 1.65o- re- 
sults in a slope ~ 1.32 ± 0.14 (shown by the dashed line in 
Fig. |4ll . The slope of this relation is consistent with our re- 
sults for black holes less massive than about IO^Mq, but for 
more massive black holes the assumption of Vopt oc V200 
breaks down. Shankar et al. (2006) use the abundance 
matching method between black holes and dark matter 
haloes to determine a relation between black hole mass and 
halo mass. They find a relation with a double power-law. At 
high halo masses (> lO^^Mo) the slope is 1.25. This agrees 
with our relation at a halo mass of M200 10^^/i-^Mq, 
but by a halo mass of M200 = 10^^/i-^Mq our relation has 
a significantly shallower slope of ~ 0.65. 

The slope of the Mbh — Af200 relation has theoretical in- 
terest, because different mechanisms for black hole growth 
are purported to predict different slopes (e.g., Wyithe & 
Loeb 2003). Our results show that the slope of this rela- 
tion is not constant, varying from ~ 1.4 at halo mass of 
Af200 ~ lO^^ft"^M0 to ~ 0.65 at halo mass of A/200 ^ 




12 15 14 

logio M200 [h-' Mq] 

Figure 4. Relation between central super-massive black hole 
mass, AfBHi and halo mass, M200- The solid line shows the re- 
lation derived assuming our relation between Vopt and V200 for 
early-types, and the A^bh — relation from Giiltekin et al. (2009). 
The shaded regions correspond to the 2a uncertainties in the re- 
lation. The solid line has a slope of ~ 0.65 at high halo masses, 
and ~ 1.4 at low halo masses, as indicated by the red long-dashed 
lines. The short-dashed line shows the A^bh — Af200 relation de- 
rived assuming V200 = Vopt = 1.65cr. This relation has a slope 
of ~ 1.32. For reference the dotted line shows a linear relation 
between ATbh and Af200- 



Table 3. Parameters of double power-law fitting formula (Eq. [Sjl 
to the y = M200 [h~^MQ] vs x = Mbu [Mq] relations in Fig.|4] 



a 


/3 


Xo 


yo 


7 


Early-types: ran 


56 logio 


Mbh = 


7.08 - J 


S.84M0 


mean 0.751 


1.633 


8.158 


12.92 


1.766 


+2o- 0.751 


1.589 


8.158 


13.19 


1.766 


-2cr 0.751 


1.678 


8.158 


12.65 


1.766 



10^^ "'h^^MQ. This might imply that different mechanisms 
for black hole growth are important at different halo masses. 
However, since the Mbh — A/star relation is roughly linear, 
the mass dependent slope of the A4"bh — M200 relation is 
driven by the mass dependent slope of the Matar — M200 
relation. This latter relation is determined by the physics 
of galaxy formation, which might not be directly governed 
by the physics of supermassive black holes. Thus the rela- 
tion between black hole mass and stellar mass may be more 
fundamental (in terms of black hole physics) than the rela- 
tion between black hole mass and halo mass. We note that 
the change in slope of the Mbh — M200 relation occurs at 
Af200 ~ lO^'^'^, corresponding to the mass scale of galaxy 
groups. 

Using a cosmological simulation for the co-evolution of 
black holes and galaxies Booth & Schaye (2010) argue that 
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black hole masses are determined by the masses of their host 
dark matter haloes. As supporting evidence for their conclu- 
sion they find good agreement between the Mbh — M200 rela- 
tion in their simulations and the observational relation from 
Bandara et al. (2009) (i.e. a slope of 1.55). Our results for 
the slope of the Mbu — M200 are consistent with the predic- 
tion of Booth & Schaye (2010) for black hole masses between 
~ lO'^ and ~ IO^Mq, but for higher black hole masses our 
results have a much shallower slope. Thus black hole masses 
cannot be universally determined by the masses of their host 
dark matter haloes via the mechanism proposed by Booth 
& Schaye (2010). 



5 COMPARISON WITH ACDM HALOES 

The majdmum circular velocity, V^nax.h, of an NFW 
(Navarro, Frenk, & White 1997) halo is reached at a ra- 
dius, rmax,h — 2.163r„2 = (2.163/c)i?200, where r_2 is the 
NFW scale radiu^, -R200 is the virial radius of the halo, and 
c = R2oo/r-2 is the halo concentration. For Milky Way mass 
haloes, c ~ 7, and thus rmax.h — 0.3-R200- The ratio between 
Knax,h and V200 is a function of the halo concentration, and 
is given by 

^ 0.465v/^7I(^, (8) 

V200 

where A{x) = ln(l + x) - x/{l + x). The factor of 0.465 is 
equal to A{x)/x with x — 2.163. 

Fig. [5] shows the relation between Vopt/V2oo and V200 
for early- (red shading) and late-type (blue shading) galax- 
ies. Note that at a fixed Vopt, uncertainty in V200 moves 
points diagonally from top left to bottom right in the 
K)pt/V2oo vs V200 plot. The dotted line shows the relation 
between Vinax,h/V2oo and V200 for ACDM dark matter haloes 
using the concentration mass relation for relaxed haloes 
(Maccio et al. 2008) in WMAP 5th year cosmology (Dunkley 
et al. 2009): logy, c = 0.830-0.098 logio(M2oo/lO^^/i"^M0). 

The green shaded region shows the variation in 
Knax,h/V2oo Corresponding to 2cr scatter (0.22 dex) in halo 
concentrations (Maccio et al. 2008). Note that the uncer- 
tainty in the zero point of the concentration mass relation 
from uncertainty in cosmological parameters (mostly ug) is 
smaller than this. 

For late-type galaxies V2.2/V200 is consistent with 
V;nax,h/V2oo, but there is a systematic trend for lower 
V2.2/V200 for low mass haloes. This trend largely disappears 
with ymax/V2oo (loug-dashed line), and indicates that the 
maximum rotation velocity measured within the optical re- 
gion of late-type galaxies is equal to the maximum circu- 
lar velocity of the dark matter halo. Such an assumption is 
commonly made in semi-analytic galaxy formation models 
(e.g. Croton et al. 2006), and when attempting to infer dark 
halo masses from optical rotation velocities (e.g. Blanton 
et al. 2008; Genel et al. 2008). Our results give this assump- 
tion empirical justification. 

For early-type galaxies the situation is more complex, 
as K)pt/V2oo is inconsistent with being constant. For mas- 
sive early-types (V200 ^ 250kms~^), Vopt/V200 decreases 

^ The NFW scale radius is the radius where the logarithmic slope 
of the density profile, dlnp/dlnr = —2, hence the notation, r_2. 




'og,o V200 [km s-'] 

Figure 5. Optical-to-virial velocity ratio vs halo virial velocity. 
The shaded regions correspond to 2a uncertainties. Early- types 
are shaded red, while late-types are shaded blue. For late-types 
the dotted line shows the relation between the maximum rotation 
velocity and the halo virial velocity. The relation between the 
maximum circular velocity of dark matter haloes and the virial 
velocity of dark matter haloes is given by the long dashed line. 
The green shaded region shows the scatter in this ratio due to 2cr 
(i.e. 0.22 dex) variation in halo concentration. 

with increasing halo mass, while for lower mass early-types 
(130 < V200 < 250kms-i), Kpt/V2oo 1.2. For halo veloc- 
ities of 1^200 — 350kms~-^, Vopt — V2oo- 

For the highest halo masses (^200 ^; 500 km s-\ A/200 ==; 
3 X lO"/i"^M0), Vopt/V200 < 1. For these halo masses the 
halo concentration c ~ 5, and thus the maximum circular 
velocity of the halo occurs at rinax,h — 300 kpc. This scale 
is an order of magnitude higher than the half light radii of 
massive early type galaxies, which are of order 10 kpc (e.g. 
Shen et al. 2003). Thus the result that T4pt/V2oo < 1 can be 
interpreted as a consequence of the galaxy half light radii 
only probing the rising part of the halo circular velocity 
curve. This mass scale also corresponds to galaxy groups, so 
it should not be a surprise that the baryons in the central 
galaxy only probe a small fraction of the virial radius. 

For lower halo masses (130 < V200 < 250 kms"^) there 
is evidence that Vopt > Vmax,h which would suggest that 
baryons have modified the potential well (either by their 
own gravity, or by modifying the structure of the dark mat- 
ter halo through adiabatic contraction). We will attempt to 
disentangle these two possibilities in a future paper (Button 
et al. in prep). However, Vopt = Vmax,h is also consistent 
with the data (within the 2a uncertainties) for this velocity 
range. 

6 COMPARISON WITH OBSERVED AND 
SIMULATED GALAXIES 

Figs. [6] & [7] shows a comparison between the V200 — Mstar, 
Vopt — Mstar and Vopt — V200 relations that we derive in this 
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Figure 6. Comparison between the virial velocity (left) and optical velocity (right) TF and FJ relations from this paper (dashed lines 
and shaded regions) and the values for observed (white symbols) and simulated (colored symbols) galaxies. The observed galaxies are 
M31 ( Klypin et al. 2002, K02, right star) and the Milky Way (Klypin et al. 2002, K02, left star; Smith et al. 2007, 307, white square, 
with 90% CL error bars; Xue et al. 2008; white diamond, with 2a error bars, for clarity Mstar has been shifted by 0.02 dex). The galaxies 
from cosmological hydro-dynamical simulations are from: Abadi et al. 2003, AOS, green pentagon; Governato et al. 2007, G07, cyan 
triangles; Naab et al. 2007, N07, red circles; Tissera et al. 2009, T09, magenta diamonds; and Piontek & Steinmetz 2009, blue squares. 
For the Vopt — Afstar relation (right panel) the shaded region shows the Icr intrinsic scatter, for the V200 — Afgtar relation (left panel) the 
shaded region shows the 2a uncertainty. Stellar masses assume h = 0.7. In the left panel the dotted lines show relations corresponding to 
integrated star formation efficiencies, egp = 10,25,50 & 100%, assuming a cosmic baryon fraction of 0.165. The MW and M31 galaxies 
are offset from both the V200 — A/star and Vopt — Afstar relations for late-type galaxies. None of the simulated galaxies fall on both the 
virial and optical TF and FJ relations. 



paper (solid lines and shaded regions) with the values ob- 
tained from cosmological simulations (solid symbols) and 
inferred for the Milky way (MW) and Andromeda (M31) 
galaxies (open symbols) . Where necessary we have converted 
the halo masses and virial velocities into the Avir = 200 def- 
inition (see Eq. [T| using the halo concentrations measured 
or assumed by the authors, or assuming Af200 = 0.83Afvir 
(where Mvir corresponds to Avir — 100) where concentra- 
tions are not specified. Table |4] contains the parameters of 
the relations between Vopt/V2oo and Vopt that we derive in 
this paper, and Table [S] contains the parameters of the data 
points in Figs. |6]&lZ] 

For the MW we adopt Vopt = 220 km s"^ (the lAU 
value) and a stellar mass of 5.2 ± 0.5 X 10^°Mq derived by 
Widrow, Pym & Dubinski (2008) using dynamical models. 
This mass is consistent with the mass of 4.8 x IO^^Mq from 
Klypin, Zhao, & Sommerville (2002) using dynamical mod- 
els, the 4.85 - 5.5 x 10^°Mq from Flynn et al. (2006) us- 
ing dynamical constraints, and 5 x IO^^Mq from Hammer 
et al. (2007) using stellar population models and a Kroupa 
IMF. 

For M31 we adopt Vopt = 260 ± 10 kms~^ based on the 
rotation curve data compiled by Widrow, Perrett, & Suyu 
(2003). We adopt a stellar mass of 9.6±0.7x 10^°Mq. This is 
consistent with the mass of 8.9 X 10"'Mq favoured by Klypin 
et al. (2002) using dynamical models, the 9.5 x lO^'^M© from 
Widrow et al. (2003) using dynamical models, and the 10.3 x 



10^" M0 from Hammer et al. (2007) using stellar population 
models and a Kroupa IMF. 



6.1 The Milky Way and Andromeda Galaxies 

Applying our results to the Milky Way (MW) (Vopt = 
220kms"^) and Andromeda (M31) (Vopt = 260kms"^) 
galaxies we predict that Vopt/V2oo = l.Htoio (2(7), with 
corresponding V200 = 198^32 kms~^ for the MW and V200 = 
235l3gkms"^ for M31. The values of Vopt/V2oo that we 
derive are significantly lower than suggested by the dy- 
namical models for the MW and M31 of Klypin, Zhao & 
Somerville (2002) which correspond to Vopt/Vjoo = 1.66 
(i.e. V200 = 137kms-i) for MW and Vopt/V2oo = 1.62 (i.e. 
V200 = 161 kms'^) for M31. 

Other measurements of the virial mass of the MW also 
imply relatively high values for Vjpt/Vjoo (i.e. low values for 
V^oo): Using high velocity stars from the RAVE survey Smith 
et al. (2007) find Mvir = lA2tol\ x lO^^M© which corre- 
sponds to V200 = 149^22 kms~^ (after converting halo defi- 
nitions) and V)pt/V2oo = 1.48lQ'2g; Using kinematics of blue 
horizontal branch stars in the MW halo, Xue et al. (2008) 
find A/vir = 1.0to;4 X lO^^M0(2o-) which corresponds to 
V200 = 132l2iknis~^ (after converting halo definitions) or 

Vopt/V200 = 

The Kahn & Woltjer (1959) timing argument estimates 
the mass of the Local Group from the age of the Universe 
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and the separation and relative radial velocity of the MW 
and M31. Li & White (2008) use a cosmological N-body 
simulation to calibrate the bias and error distribution of 
the Timing Argument estimators of the masses of the Local 
Group. They derive a combined virial mass of the MW and 



0.4 



+4.98 



X 10'^ (90% CI). 



M31 of M200 = 5.25 

To interpret this mass we need to divide it between the 
MW and M31. We consider two simple ways to do this: Case 
1: Equal halo masses. M200 = 2.64li;72 x lO^^M©, or V200 = 



2.64Ii?^ X 10" 

199.11*^:^ kms-\ which implies Vopt/V200 = l-Wto.tl for 



the MW, and V;>pt/V^2oo = 1.3lj 



for M31. Case 2: Equal 



K)pt/V2oo. This implies a halo mass ratio of (260/220) = 
1.65, and results in Kpt/V2oo = l-'iltoH (90% CI) for both 
the MW and M31. Thus the timing argument gives Vopt / V200 
in better agreement with our determination, but given the 
large measurement uncertainties it is also consistent with 
K>pt/F2oo ^ 1.6 for the MW. 

If the discrepancy between the l/opt/V2oo ratios as de- 
termined by our analysis and that measured for the MW 
and M31 holds under further study, it would imply that the 
MW and M31 do not live in typical dark matter haloes for 
their optical rotation velocity. This would also imply that 
there is substantial scatter in the V'opt/V2oo ratio. Scatter in 
Kjpt 7^^200 is expected. For example, the analytic models of 
Button et al. (2007), which are calibrated against the TF 
and size-luminosity relations predict a la scatter of ~ 0.05 
dex in this ratio. However, reconciling yopt/V2oo — 1.6 with 
our measurement of K)pt/V2oo — 1.1 requires a high sigma 
outlier. If we require the MW and M31 to be only slightly 
atypical, a larger scatter would be needed, of at least 0.1 dex. 
However, a large scatter in K)pt/V2oo may be difficult to rec- 
oncile with the small scatter (~ 0.05 dex in velocity) in the 
observed {Vopt — Mstar) TF relation (Courteau et al. 2007; 
Pizagno et al. 2007). 

If the MW and M31 have atypical V'opt/V200, the left 
panel of Fig. [S] shows that the MW and M31 also have 
atypical stellar masses for their halo velocities, with the 
stellar masses being higher than average. In terms of the 
Vopt - Mstar TF relation, the MW and M31 are also atypi- 
cal, being offset to high velocities by more than la (in terms 
of the intrinsic scatter). The fact that the MW does not 
fall on the TF relation (including the /-band luminosity, 
stellar mass, and baryonic mass variants) has been noted 
previously (Flynn et al. 2006; Hammer et al. 2007). Ham- 
mer et al. (2007) found that M31 does fall on the TF rela- 
tion. This apparent conflict with our result is due to the 
different velocity definitions used. Hammer et al. (2007) 
adopt Vfiat = 226kms~^ for M31, whereas we adopt V2.2 = 
260kms~^ (based on the rotation curve data compiled by 
Widrow et al. 2003). 
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Figure 7. Comparison between the observed Vopt/Vioo - Vopt 
relation from this paper (shaded regions correspond to 2cr un- 
certainties) and the values of observed (open symbols) and sim- 
ulated (solid symbols) galaxies. The observed galaxies are M31 
(Klypin et al. 2002, K02, right star) and the Milky Way (Klypin 
et al. 2002, K02, left star; Smith et al. 2007, S07, white square, 
with 90% CL error bars, for clarity Vopt is shifted by -O.Oldex; Xue 
et al. 2008; white diamond, with 2cr error bars). The galaxies from 
cosmological hydro-dynamical simulations are: Abadi et al. 2003, 
AOS, green pentagon; Governato et al. 2007, G07, cyan triangles; 
Naab et al. 2007, N07, red circles; Tissera et al. 2009, T09, ma- 
genta diamonds; and Piontek & Steinmetz 2009, blue squares. See 
text for further details. 



Table 4. Parameters of double power-law fitting formula (Eq. [Sjl 
to the y = Vopt /V200 vs x = Vopt relations in Fig. [7] 



xo 



7 



Early-types: rangclogj^Q Vopt = 169 — 470 km s ^ 
mean 0.009 -1.156 316.4 1.078 6.993 
+2ct 0.009 -1.214 316.4 1.326 6.993 
-2a 0.009 -1.098 316.4 0.876 6.993 

Late-types: range logj^Q Vopt = 91 — 265kms~^ 
mean 0.407 -0.186 180.1 1.099 3.559 
+2cr 0.348 -0.186 180.1 1.321 3.559 
-2cr 0.585 -0.186 180.1 0.886 3.559 



6.2 Galaxies from Cosmological Simulations 

In Figs, [g] & [3 we also show the V2OO — Mstar, Vopt — Matar 

and Vopt — V200 relations of galaxies formed in hydrody- 
namical cosmological simulations: Abadi et al. (2003, green 
pentagon); Governato et al. (2007, cyan triangles); Naab 
et al. (2007, red circles); Tissera et al. (2009, magenta dia- 
monds); and Piontek & Steinmetz (2009; blue squares). For 
the simulations of "disk" galaxies of Abadi et al. (2003), and 
Governato et al. (2007), we use Vopt measured at 2.2 disk 
scale lengths (/-band for Abadi et al. 2003, K-hand for Gov- 



ernato et al. 2007). These simulations predict V)pt/V2oo — 2, 
which is even more discrepant with our results than the es- 
timates for the MW and M31. However, as shown by But- 
ton & Courteau (2008) the simulated galaxies from Gov- 
ernato et al. (2007) do not fall on the TF relation (V2.2 
vs /-band luminosity), being offset to high rotation veloci- 
ties. This offset is also apparent in the Vipt — A/atar relation 
as shown in Fig. The cause of this offset is not clear, 
but it is likely a combination of insufficient numerical res- 
olution which leads to artificial angular momentum losses 
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Table 5. Velocity and mass parameters of the observed and simulated galaxies shown shown in Figs.|6]&[7] 



Galaxy ^^200 


V200 




Vopt 


Vopt/V200 




Mstar 






[km/s] 




[km/s] 








[lO^OMnl 




This paper 


















MW 18.0(+14.9, -7.5) 


198(+44, 


-32) 


220 


1.11(- 




-0.20) 


5.2(+0.5, 


-0.5) 


M31 30.1{+23.6, -12.4) 


235(+50, 


-38) 


260 


1.11(- 


h0.21 


-0.19) 


9.6(+0.7, 


-0.7) 


Observed: timing argument, Li & 


r. White (2008) 














MW+M31 52.5(+49.8, -33.0) 


















Observed: halo stars, Xue et al. (2008) 
















MW 8.16{+4.90, -3.26) 


132(+22, 


-21) 


220 




h0.31, 


-0.24) 


- 




Observed: high velocity stars. Smith et al. (2007) 














MW 11.8(+9.4, -4.5) 


149(+32, 


-22) 


220 


. 


hO.25, 


—0.26) 






Observed: dynamical models, Klypin et al. (2002) 














MW 8.59 


137.2 




228 


1.66 






4.8 




M31 13.7 


160.5 




260 


l.DZ 






8.9 




Simulated: Abadi et al. (2003) 


















5.6 


134 




260 


1.94 






10.3 




Simulated: Governato et al. (2007) 
















DWFl 1.38 


75 




165 


2.21 






1.38 




MWl 10.0 


144 




280 


1.94 






4.60 




GALl 26.7 


200 




330 


1.65 






13.6 




Simulated: Naab et al. (2007) 


















Al 16.9 


167.8 




270 


l.Di 






21.7 




A2 17.9 


171.0 




232 


l.oO 






22.8 




CI 12.5 


151.7 




258 


1. / U 






17.2 




El 13.4 


155.3 




307 


1.98 






19.0 




Simulated: Tissera et al. (2009) 


















Aq-A-5 11.0 


167.9 




251 


l.OU 






8.11 




Aq-B-5 5.2 


130.8 




166 


1 0'7 
i.Z / 






3.47 




Aq-C-5 11.8 


171.9 




257 


l.OU 






8.12 




Aq-D-5 10.9 


167.4 




229 


1 'i'7 






6.04 




Aq-E-5 7.9 


150.3 




224 


1 /1 






6.81 




Aq-F-5 6.7 


142.3 




195 


i.o 1 






7.38 




Simulated: Piontek & Steinmetz 


(2009) 
















DM-hrl 0.82 


70.7 




105.4 


1 /in 






0.64 




DM-hr2 1.53 


87.0 




136.9 


1.57 






2.11 




DM-hr3 2.18 


97.8 




156.1 


1.60 






2.36 




DM-hr4 3.00 


108.9 




186.5 


1.71 






3.61 




DM-hr5 3.61 


115.8 




199.0 


1.72 






4.16 




DM-hr6 4.28 


122.5 




194.8 


1.59 






4.48 




MW-hr 7.48 


147.6 




207.5 


1.41 






6.41 




DM-mr7 14.4 


183.5 




319.6 


1.74 






16.7 





(Kaufmann et al. 2007), and/or insufficient feedback which 
results in baryon fractions that are too high (Button & van 
den Bosch 2009), and/or too much adiabatic contraction of 
the haloes. 

The more recent simulations of disk galaxies from Pi- 
ontek & Steinmetz (2009) predict values of V2. 2/1^200 — 
1.6 ± 0.1. Although these are lower than the V2.2/V200 from 
Abadi et al. (2003) and Governato et al. (2007), they are 
still highly inconsistent with our measurements. Fig.|6]shows 
that while the simulations from Piontek & Steinmetz (2009) 
fall on the Vopt — Afstar relation at low masses, the simulated 
galaxies have too much stellar mass at fixed halo velocity (or 
halo mass), especially at low masses. 

This demonstrates that while cosmological simulations 
of disk galaxies have made great progress in producing galax- 



ies that fall on the V200 — A^star and Vopt — Mstar relations, 
they still have been unable to produce galaxies that simul- 
taneously reproduce both of these relations. Reproducing 
these relations will provide a key test for galaxy formation 
models. 

The red circles show the simulations of Naab 
et al. (2007), with parameters taken from Johansson 
et al. (2009). The simulations produced early-type galaxies 
(spheroids with no disk component). Three galaxies have 
Vopt/V2oo > 1.6, which is higher (by more than 2a) than 
our derived value. Galaxy A was re-simulated with 8 times 
more particles (200'^) resulting in Vopt = 232kms~^ and 
V)pt/V2oo = 1.36, this galaxy is shown with a larger red cir- 
cle. This simulation has V]pt/V2oo within la of our results. 
This shows that numerical resolution is still an important 
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issue for cosmological simulations that wish to resolve the 
internal structure of galaxies. While the simulations from 
Naab et al. (2007) can produce a relatively low value of 
Vopt/V2oo, Fig. E] shows that this gala:xies does not fall on 
the V200 — Mstar or Vopt — Mstar relations for early types (red 
shaded regions). 

The magenta diamonds show the simulations of Tissera 
et al. (2009), where we have estimated Vopt from their Fig- 
ure 9. These simulations have K)pt/V200 — 1.4 ±0.1, which is 
slightly higher, but consistent within the uncertainties with 
our results for early-type galaxies. The simulations of Tis- 
sera et al. (2009) have only 10^ particles per galaxy, and 
thus (based on the resolution tests of Naab et al. 2007) the 
Vopt /V200 may be biased high by insufficient numerical reso- 
lution. Fig. [6] shows that the Tissera et al. (2009) simulated 
galajcies also have stellar masses that are too high for their 
halo and optical circular velocities. 

The fact that current cosmological simulations produce 
galaxies with too many stars has also been shown by Guo 
et al. (2010). This should result in an increase in V^pt/V2oo, 
and thus it is plausible that the high values of T4pt/V2oo 
found in the simulations are just the result of the stellar 
masses being too high. Ifowever, having the correct amount 
of stars does not guarantee that a model galaxy has the cor- 
rect 14pt/V2oo. For example, the two most massive galaxies 
from Governato et al. (2007) have roughly the correct stel- 
lar masses for their halo masses (for late- types), but they 
have K)pt/V2oo — 2 which is much too high. For a fixed 
V200 and Mstar the 14pt/V2oo ratio also depends on the size 
of the galaxy, the halo concentration and on the response 
of the halo to galaxy formation (e.g. Mo, Mao, & White 
1998; Button & van den Bosch 2009). Smaller galaxies, 
higher halo concentrations and halo contraction all result 
in higher Vopt/V2oo. Thus in order to determine the origin 
of the V^pt/1^200 ratios, the sizes of galaxies are an essential 
observational constraint. 



7 SUMMARY 

We combine measurements of halo virial masses from weak 
lensing, satellite kinematics, and halo abundance match- 
ing with the TuUy-Fisher (1977) and Faber-Jackson (1976) 
relations to place constraints on the average relation be- 
tween the optical (Vopt) and virial {V200) circular veloc- 
ities of early- (Vopt = V50 = 1.65(T(i?5o)) and late-type 
(Vopt = V2.2 = Vot(2.2i?d)) galaxies at redshift z ~ 0. We 
summarize our results as follows; 

• The stellar mass to halo virial mass fractions of late- 
type (blue/disk dominated) galaxies increase from ~ 0.015 
at a stellar mass of A'/star = lO^ */i"^M0, to ~ 0.043 at a 
stellar mass of Matar ^ lO" °/i"^M0 (assumiug h = 0.7). 
These correspond to integrated star formation efficiencies, 
esF = Mstar/(/barM2oo), of ESF 9% aud esF ^ 26%, 
respectively (assuming a cosmic baryon fraction, /bar ~ 
0.165). After accounting for cold gas, the galaxy formation 
efficiencies, ecF = (Mstar + Mgas)/(/barAf2oo), are eoF =i 
17% and eoF — 33%, respectively. 

• The central galaxy formation efficiencies of early-type 
(red/bulge dominated) galaxies reach a peak of ~ 12% at 
Mstar = W^°-^h~^MQ. The efficiency drops to ecF ^ 2.8% 
at Mstar = 10"-*/i-''Mq. 



• For late-type galaxies V2.2/V200 is roughly constant, and 
close to unity, over the range of stellar masses or rotation 
velocities that we probe: (10^'^^ < Mstar < lO"/i"^M0, or 
90 < V2.2 < 260kms'i). 

• For late-type galaxies, the maximum circular velocity 
of ACDM haloes, Vmax.h, is consistent with being equal to 
the maximum observed rotation velocity, Vmax. This result 
is not reproduced by current cosmological simulations of 
disk galaxy formation (e.g. Governato et al. 2007; Piontek 
& Steinmetz 2009). 

• Cosmological simulations have been unable to form 
galaxies that fall on both the V200 — Mstar and Vopt — Mstar 
relations. These relations provide a strong and simple test 
for galaxy formation models. 

• For early- type galaxies V50/V200 is not a constant, 
which rules out most early-type galaxies from being glob- 
ally isothermal (i.e. constant circular velocity from the op- 
tical half light radius to the virial radius of the halo) . 

• Early-type galaxies have V50/V200 — 1 only for Vopt — 
350 kms"'^ (i.e. a ~ 215kms~^). This is consistent with the 
strong plus weak lensing result of Gavazzi et al. (2007) which 
inferred that massive early-types with a ~ 240 km s^^ have 
close to isothermal global density profiles. 

• There is some evidence that for early-type galaxies in 
lower mass haloes (V200 250kms~^) the V50/V200 ratio 
is higher than the Vmax,h/V2oo ratio, which would indicate 
that baryons must have modified the potential well, either 
through their own gravity, or by adiabatic contraction of the 
halo. 

• For early-type galaxies in higher mass haloes (V200 ^ 
500 km s"'^), V50 /V200 < 1, indicating that the half light radii 
of the galaxies are much smaller than the NFW scale radius. 

• The mass dependence of the V50/V200 ratio for early- 
types implies that the slope of the black hole mass - halo 
mass relation varies with mass. Specifically it varies from 
~ 0.65 at high halo masses (~ lO^^-^M©), to ~ 1.4 at low 
halo masses (~ W^^Mq). The shallow slope at high mass is 
at odds with that predicted by the self-regulating feedback 
model of Wyithe & Loeb (2003) and the cosmological simu- 
lation of Booth & Schaye (2010) which both predict a slope 
of ~ 1.5. 

• For the Milky Way (MW) (Vopt = 220 kms"^) and An- 
dromeda (M31) (Vopt ~ 260kms~^) galaxies our results pre- 
dict that 



Vopt/V2oo = 1.11 (+0.22, -0.20, 2a) (9) 



This value is significantly lower than those suggested by the 
dynamical models for the MW and M31 of Klypin, Zhao & 
Somerville (2002) which correspond to Vipt/Vaoo = 1.66 for 
MW and V)pt/V2oo ~ 1-62 for M31. Other measurements of 
the virial mass of the MW using kinematics of halo stars 
also infer relatively high values: Vopt/V2oo = 1.48lo (90% 
C.I.) (Smith et al. 2007); Vopt/V2oo = l-67l°ji (2o-) (Xue 
et al. 2008). If both our relation and the dynamical mod- 
els are correct this discrepancy would imply that the MW 
and M31 do not live in typical dark matter haloes for their 
optical rotation velocity. 
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